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Abstract. We propose an extension of matrix element plus parton shower merging at tree level to preserve 
inclusive cross sections obtained from the merged and showered sample. Implementing this constraint gen- 
erates approximate next-to-leading order (NLO) contributions similar to the LoopSim approach. We then 
show how full NLO, or in principle even higher order, corrections can be added consistently, including con- 
straints on inclusive cross sections to account for yet missing parton shower accuracy at higher logarithmic 
order. We also show how NLO accuracy below the merging scale can be obtained. 



' PACS. 12.38.Bx Perturbative QCD calculations - 12.38.Cy Summation of QCD perturbation theory 
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(D 1 Introduction 

fi: 

, Parton shower event generators, [D{3], to name only few of the most recent general- purpose simulations, are by now 
j^. ■ indispensable workhorses of experimental as well as theoretical studies for comparing (standard model) theoretical 
(■ — 1 predictions to measured observables in a most detailed way. Within recent years, tremendous progress has been made 
\Q . in improving the approximations underlying those simulations by exact calculations in perturbation theory. Starting 

■ from simple matrix element corrections on the one hand, improving the hardest parton shower emission to be driven 
I/"") | by the exact tree level matrix element [H[S] , subsequently tree level matrix elements for n > 1 hardest emissions have 

■ been included consistently [6lll3]. These approaches of merging multileg matrix elements at tree level, i.e. leading 
t—{ order (LO), have been accompanied by efforts of matching parton shower simulations and perturbative calculations 
£N| ■ at next-to-leading order (NLO), including both matrix elements for additional hard emission at tree level as well as 
t~ I ■ virtual, one-loop corrections to the hard process of interest, [T4Ul9j . On the other hand, parton shower algorithms 

\ themselves have evolved from crude approximations of multiple parton emission to more and more refined and precise 
• i-H ■ tools [20H26] , eventually aiding attempts of combining fixed-order calculations with shower resummation in a most 
1 consistent way. 

NLO matching has been used to improve tree level merging algorithms, [2TH50] . but only recently, taking full 
advantage of developments in both merging, matching and shower algorithms, first experience has been gained in 
combining NLO QCD corrections to multijet final states of different multiplicity with subsequent parton shower 
emissions [3"TM34| . The aim of such approaches is to obtain an event sample which will provide LO (NLO) precision for 
observables controlled by hard n-parton emission as long as exact QCD corrections are included for up to n partons 
at tree level (n partons at one loop and n + 1 partons at tree level), while including resummation by means of shower 
simulation at whatever accuracy is provided by these algorithms. 

While merging approaches at tree level have been proven to be amazingly powerful in describing experimental 
data, and can now be considered to be well understood at a theoretical level (particularly with respect to logarith- 
mic dependence on the merging scale separating hard, matrix element driven emissions from softer parton shower 
emissions), their generalization to include NLO corrections yet suffers from impact of the merging scale at a level of 
logarithmic accuracy which by no means can be provided by existing shower algorithms. The most striking signs of 
this dependence are expected in inclusive cross sections as predicted by such simulations, and proposals to cure the 
problem by explicit input of resummation results at higher accuracy have been made |35j . In the present contribution 
we try to approach this problem from a rather pragmatic point of view by setting up a formalism to systematically 
include higher oder corrections in parton shower simulations, while satisfying constraints to obtain the proper inclusive 
cross sections at the respective order of perturbative calculations available to the merging algorithm. 
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2 Setting the scene 

We will consider a generic parton shower with an evolution variable q (all other kinematic variables of a splitting are 
suppressed for the sake of readability), which is driven by splitting kernels P((j) n ,q). Here and in the following, 4> n 
denotes a partonic configuration (phase space point) containing n additional partons with respect to the lowest order 
process of interest, and P(4> n ,q) determines the dynamics of emission at a scale q off the partonic system <j> n . The 
Sudakov form factor associated to P(4> n , q), evolving from a hard scale Q to a soft scale q, is given by 



A n {q\Q) = exp 



Q 



dk 



l>n+l 

b n dk 



P(4>n,k) 



(1) 



with the phase space Jacobian d^„+i / d<j) n depending implicitly on the shower kinematic variables. Partonic config- 
urations <j) n are determined according to differential cross sections da(<p n , q n \ ■ ■ ■ \qo), where the sequence of scales 
q-n < <Zn-ij ■■■,9o is either directly determined from shower evolution, or is assigned by a clustering procedure corre- 
sponding to the inverse of a possible shower evolution 4>q — > <f> n , if the cross section is determined by exact matrix 
elements!]] In the case of exact matrix elements, reweighting to account for couplings and PDF factors evaluated at the 
clustering scales may as well be included in this notation. For the sake of readability we will skip the whole sequence 
of scales and denote da(<p n ,q n \...\qo) — da((j) n ,q n ). The parton shower action on cross sections for <j> n events is given 
bj0 



PS M [da((f> n ,q n )] = da(<f> n ,q n )A n (iJ,\q n ) + PS A 



da(4> n ,q n )- 



P^,{4> n ,q n+ i)A n (q n+ x\q n ) 



(2) 



The first contribution are events where no further radiation has been generated (/x denotes the parton shower cutoff 
and the splitting kernels P /i (^!> n , q n +i) — Pi&n, qn+\)d{q n +i — p) vanish for scales below the cutoff), while the second 
contribution corresponds to at least one emission. 

Considering parton showering off the lowest order tree level cross section, dcr^ \(j)Q, qo), we find that the first N 
iterations of the parton shower action give rise to 



N-l 



PS, 



k=0 



da^((f> , q ) = da^(cj> ,q ) £ -^P^fa^, q k ) ■ ■ • F M (0 O , qi )A k ^\q k \ ■ ■ ■ \q Q ) + 



PS„ 



da^(cf>o,q 



d<j) 



N 



Pu 



N 



-i,gjv) • • -Pn((f>o, qi)A N -i(q N \ ■ ■ ■ \q ) 



, (3) 



where Af.(/J,\qk\ ■ • ■ \qo) = Ak(p\qk) • ■ • Ao(q%\qo). Considering the exclusive and inclusive cross sections for producing 
exactly n or at least n partons, respectively, we have 



> n 



da {0) ((/) ,q )'^j } -P t _ t (<j) n ^ 1 ,q n ) ■ ■ ■ P^o, q 1 )A n (fj,\q n \ ■ ■ ■ \q ) 
acpo 

d<7 {0) ((t) ,q Q )'^p-P t _ l (<j) n ^i,q n ) ■ ••P /1 (<^o,gi)Ai-i(gn| • • • \Qo) 
acpo 



Note that the latter expression is a direct implication of the parton shower being unitary. In particular, the total 
inclusive cross section is not altered by the parton shower and determined by the lowest order 'input' cross section 
du' ' , as would be expected from the very definition of an inclusive cross section. 

This property is a direct result of the fact that the emission contribution of the parton shower at scale q is an exact 
differential of the no emission probability down to q, 



qk-i 



-P(4>k-i,qk)A k -i(q k \qk-i) = 1 - A k -i(q\q k -i) 



(4) 



-i,q k )=P{ 



d0 fe _idgfc 

k -i,q k )9(p — qk)6(qk — A*) with p > fi, and associating PS M | p as the associated action 



Introducing P J[1 |p( 
we also have 

PS„[-] = PS M | P [PS P [-]] , (5) 

i.e. we can always split the shower evolution into two (or more) scale ranges. We illustrate the main formulae discussed 
in this section in a diagramatic manner in fig. [T] 

Note that ordered histories may not always exist. In these cases we can assume that, for example, the history with the 
longest ordered history from smaller to larger scales has been chosen. This will not pose a problem for the formalism outlined 
here, provided we set Sudakov form factors with unordered scales equal to one, A n (Q\q) — 1 for Q > q. Also note that a 
clustering down to the lowest order process may not always be possible, in case of which we terminate the clustering sequence 
at the last possible clustering step. 
2 For the following note that PS[-] is linear. 
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(a) (b) (c) 

Fig. 1. Main properties of parton showers illustrated in a diagrammatic way. Each box with a black frame denotes a n-parton 
contribution, where n runs from top to bottom. Rows need to be summed to obtain the total cross section driving a given 
parton multiplicity. The parton shower evolves along the vertical red lines from larger (top) to smaller (bottom) scales. Double 
red lines denote the Sudakov factors associated to a given evolution interval, horizontal arrows branching off denote emission 
at a certain scale. With the parton shower cutoff at the bottom of each box, part (a) shows exclusive cross sections for zero, 
one and two partons emitted, cf. eq. [3] Part (b) illustrates the shower cross section integrated over the contributions of more 
than two emissions, making use of eq. [4] the form of inclusive cross sections is obvious, as each column of two boxes integrates 
to zero. Part (c) illustrates how the shower splits up into two evolution intervals, cf. eq. [5] 



3 Tree-level merging 



Having analyzed the exclusive and inclusive cross sections for n-parton production by the shower, we will now review 
how these predictions can be improved by including higher multiplicity tree-level matrix elements. Clearly, the desired 
result in terms of exclusive n-parton states is to replace the product of splitting kernels approximating an exact tree- 
level matrix element by the full answer. To this extent, in the presence of a maximum number of N additional partons 
to be described by matrix elements, we are therefore to satisfy a merging condition 



N-l 

d^:; SCd ] =J2 d ^ ] (<l>k,qk)MM---\qo) +PS P [d*W(<t> N ,q N )A N ^(q N \. ■ .\q )\ . (6) 



fe=0 



As is the case for plain parton shower emissions, there are no emissions with a scale below the shower cutoff /i, which 
we here apply to regulate the divergences in the additional tree-level matrix elements considered. For this reason we 
have introduced dcr^ O) (0„, q n ) = dcr^(0 n , q n )6(q n - pi). 

It is straightforward to show that the merging condition can be satisfied by 

da™ d = d<j(^,q ) +E U^H^qk) - ^-J%(&_x, q k )da^ (cj> k ^, q^)) A k ^{q k \ ■ ■ ■ \q ) . (7) 

k=l ^ 9k I J 
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Note that we are free to introduce a merging scale p > p, by acting the parton shower on da^ c p scd . The result is similar 



to eq. IH upon replacing 



da^(<f> k ,q k ) 



fc) Qk) 



1=1 



k-i 



(8) 



i.e. emissions at scales smaller than p are purely driven by shower dynamics, whereas above the exact tree level matrix 
element is used. The merging scale at this point can be viewed as just an efficiency tweak for not having to evaluate 
exact matrix elements in a region where they are well approximated by the shower. 

For later purposes, and to make the connection to CKKW-type merging [6 f 3 explicit, we rewrite the shower cross 
section in presence of the merging scale p as 



N-l 



PS L 



da 



merged 
N,p 



X>S M | P da^(4> k ,q k )A k (p\q k \---\q ) 



d<rf( 



*>N, qN)A N ^i(q N \ 



(9) 



fc=0 



Note that, except for the highest matrix element multiplicity N present, shower emissions are confined to happen 
below the merging scale. Using eq. [5j the previous result is simply given by 



PS L 



da 



N,p 



PS 



'JV-l 

E 

fe=0 



da p \4> k ,q k )A k (p\q k \---\q )+PS p \dtr<®(<t> N ,q N )A N - 1 (q If \---\q ) 



(10) 



In the merged sample, cross sections for exclusive n < N parton configurations above the merging scale p are determined 
by the respective tree level matrix elements including the proper Sudakov suppression, while exclusive cross sections for 
n partons down to the shower cutoff are determined by a mixture of matrix elements and shower splitting functions 
depending on which part of the relevant scale sequence is below or above the merging scale. The proper Sudakov 
suppression is as well retained in the latter case. 

Note that the merged cross section and its showered counterpart, eq. 1101 cover several approaches to tree level 
merging so long as the merging resolution coincides with the parton shower resolution in case of which no issues with 
truncated/ vetoed showers [HJH2] do appear. Throughout the paper we will assume that this is always the case; a 
mismatch between merging and shower resolution does not pose a conceptual problem to the algorithm proposed in 
the following. In the presence of a merging scale p, the formalism outlined here covers the standard merging algorithms 
of reweighted tree level matrix elements with showering below the merging scale. In this case, upon putting P p — > 
except for emissions off the highest multiplicity, the shower subtractions are absent in dcr mergedjv '' > . If no merging scale 
is present, iterated matrix element corrections are contained in the master formula upon replacing the shower splitting 
kernels for emissions off a system of up to N — 1 partons by the respective ratio of tree level matrix elements. In this 
case, the merged cross section coincides with the lowest order cross section. To work in a most generic setup, we will 
stay with the general solution to the merging condition, which does not require any modification to the shower acting 
downstream the hard processes obtained from the merged cross section. 

For convenience, let us introduce the functional 



PS 1 [da p ((/) n ,q n )} = 



do>(0n, q n ) _ d^n+l 
A n (p\q n ) d(/) n 



u<7n+l) 



do>(0n, gn) 
A n (p\q n+1 ) 



(11) 



which satisfies 



PS^ [PS M 1 [dcr M (0„,g n )]] = da p ((f) ni q n ) and PS M [PS p 1 [dcr p (0„, q n )}] = PS^ [da p (<f> n , q n )] . (12) 
Using this definition, we can rewrite the merged cross section as 

AT-l 

^PS; 1 \da p °\<t> k ,q k )A k (p\q k \---\q )\ + da p °H^ N ,q N )A N ^( qN \ ■ • • \q ) , (13) 



da™ p gcd - 



k=0 



which summarizes the solution to the tree level merging condition in a most transparent way. Being of purely formal 
use for the present letter, it is worth noting that PS" 1 facilitates the subtraction of parton shower contributions 
above the scale p at all orders. Turning this observation around, it can actually be used as a generating functional for 
subtractions needed to satisfy fixed-order matching conditions upon expanding to the respective order. Indeed, 



PS/ [dc^ (</>„,<?„)] 



do-,, 



,q n ) [ 1+ I dq n+r 



b n dq 



P, 



71+ 1 



D Qn+l) 



-Pu 



n i Qn 

) + 0(a 2 s ) (14) 




Fig. 2. The left panel shows the solution to the merging condition for N — 2, eq. 1101 with inconsistencies in inclusive cross 
sections made explicit. The right panel shows the cross sections obtained from the merged sample upon including counter terms 
to restore inclusive cross sections as given by eq. 1151 Here, emissions calculated according to exact tree level matrix elements 
are made explicit by arrows branching off blue lines. 



is readily identified as the corrections required for NLO matching. Note also that the action of PS" 1 can actually be 
implemented in a Monte Carlo simulation by performing the veto algorithm from smaller to larger scales while using 
a simple prescription to exponentiate —P((ft n ,q) at the expense of introducing weighted events |36| . 



4 Inclusive cross sections 



Let us now turn to inclusive cross sections. We are free to split up the phase space integration into a region of 
scales below the merging scale p, and one above. Performing the integration over smaller scales first (or, equivalently, 
considering jet cross sections in a jet measure exactly corresponding to inverting the parton shower), we effectively 
remove the shower action PS M | p from eq. 1101 

Looking at inclusive cross sections of at least n < N jets, we face essentially a non-unitary evolution due to the 
matrix elements not coinciding with products of splitting kernels anymore. While for the highest multiplicity we find 
the same pattern as present in the parton shower by the very definition of the merging condition, 

>N da ( f\<j) N ,q N )A N „ 1 (q N \---\q ) , 

the result for one parton less does not anymore resemble the functional dependence present in the shower 3 

> N - 1 da(°\4>N-i,qN-i)AN-2(qN-i\ ■ ■ ■ |<7o) + 

p- 1 ( d<rj?\<l> N ,q N ) d^v , , (ow, ,\ A v 

/ d 1 N j "77 a — PpVPN-uWatr), 'WN-i,qN-i) A N _i(q N \ ■ ■ ■ \q ) . 

J p V aq N d(p N -idq N p J 

In order to restore the form of the cross section obtained by the shower, we need to subtract the second term 
obtained for the inclusive N — 1 jet cross section. This pattern then continues, and one possibility to solve the issues 



3 We have added and subtracted the contribution from a shower emission off the N — I parton state, where the positive term 
has been absorbed by the integration over the scale gjv to yield the Sudakov form factor in the first line. 
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with inclusive cross sections is then given by the replacement 



da 



merged 
N,p 



j merged 
d(T N,p 




k^k+ijda^^k.qk) A k (q, 



fe+i 



W) 



A k (q k+l \q k ) A k -i(q k \ ■ ■ ■ \q ) 



+ da^(4> N ,q N )A N ^(q N \---\q Q ) . (15) 

The solution to the merging condition as discussed in section [31 as well as the results obtained from the merged sample 
including corrections for inclusive cross sections are depicted in hg. [2l One remark is in order here: while we have now 
achieved to retain the expected form of inclusive cross sections (including the constraint that the total inclusive cross 
section is given by the lowest order input cross section) the exclusive cross sections with n < N partons above the 
merging scale take a seemingly different form, 



PS 



• Qn) 



<2n 



dq„ 



da 



(0), 



t>n+l, Qn+l) 



dq 



n+l 



n I ' ' ' I Qo ) 



It is clear that such a change will occur when trying to restore inclusive cross sections. We stress however that the 
form encountered here does not pose a problem in terms of logarithmic accuracy, as long as the shower kernels are 
truly capable of reproducing the singly-unresolved limits of the tree level matrix elements. In this case the dominant 
contributions of integral in the second term will render those singularities to precisely integrate to da^\(f> n , q n )(l — 
A n (p\q n )), restoring the expected Sudakov suppression. 



5 Injecting NLO corrections 



Having improved the parton shower by tree level matrix elements, including constraints on inclusive cross sections, 
we can now turn to the case of including NLO corrections, which we assume are here available up to M additional 
partons with respect to the lowest order process. 

Before discussing this in greater detail, let us consider the respective contributions obtained from the merged 
sample with corrections for inclusive cross sections. The cross sections to be analyzed are the contributions of exactly 
n partons above the merging scale (integrated over phase space below p), which we will have to correct to be driven 
by the respective NLO cross sections. The exclusive cross sections for n partons above the merging scale are given by 



j nLO a 

ao ri,cxcl — ^r. 



n 5 Qn ) 



Qn 



dq n 



da 



(0), 



i+lj Qn+lj 



dg„ H 



A n (q n+1 \q n ) 



(16) 



Note that this does contain the Sudakov form factor for obtaining at least n partons above the merging scale. By the 
discussion in the previous section, the contribution in brackets will also provide for the remaining Sudakov factor down 
to the merging scale. Note also, that upon expanding the contribution in brackets to first order in a Sl we obtain the 
approximated, exclusive NLO cross section for n jets above a resolution p as given by the LoopSim prescription [37], 
which motivated the label nLO. 

This observation directly implies that, if we add a correction cross section given by the inclusive NLO corrections, 



ps; 



daW 



ii Qn 



^ d<r(°)(0 n+1 , gn+1 ) fl , 
dq n + 1 ( 



p) I A n -i(q n 



dg„+i 

we will obtain NLO accuracy for exclusive n parton configurations above the merging scale 

dcr (o) (0 I1+ i,g„ + i) 



da 



NLO 



190 ) 



) do-W q n ) + da^(<p r , 



dq rl 



dq„H 



6(q n -p) + 0(ai) 



(17) 



The Sudakov weight attached to the correction cross section is not running down to the merging scale, as one would 
naively expect. In this case, a double counting of logarithms of the merging scale would have happened between the nLO 
cross section and the full one-loop correction. Indeed, the role of the correction we apply here is to precisely replace 
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the first order virtual approximation stemming from expanding the Sudakov form factor by its exact counterpart. 
Fulfilling this condition is at the heart of the merging efforts at NLO presented so far, [3lTl34| . 

Let us stress the importance of maintaining inclusive cross sections in the case of NLO corrections. The problematic 
contribution by which the inclusive cross section differs from its expected value in the case of tree level merging, as 
discussed in section. [H is given by 



S (d(T n ,i nc l) 



da 



(0), 



n+l 



-1, Qn+l) 



dq n ^ 



d<j) n dq n+1 p 



^qn+^da^^n^qn) A n (q 



n+l | 



l?o) 



This cross section could well be expected not to contribute logarithmically enhanced terms, provided the shower is a 
good approximation to singly unresolved limits of tree level matrix elements. Upon replacing the LO exclusive cross 
sections above the merging scale by their NLO counterparts, the above argument does not apply anymore owing to 
the case that NLO corrections to the shower splitting kernels are not considered. We therefore expect violations of 
inclusive cross sections at a level of logarithmic approximation which is not covered by the shower anymore. Note 
that if corrections are only available to the lowest order process n = 0, then the constraints on merged inclusive cross 
sections will in this case ensure that also the inclusive NLO cross section is preserved. This is however not the case 
starting from adding NLO corrections to the n = 1 processes. 

As for the tree level case, we can however provide correction terms similar to the LoopSim nNLO corrections. In 
the presence of tree level matrix elements for up to N additional partons, and one-loop corrections for up to M < N 
additional partons we then arrive at the merged cross section, 



JV-l 



da 



merged 
N,M,p 



E *V 



fe=0 



daW(<j> k ,q k ) 



dq k -\ 



da 



(0), 



i>k+i, qk+i) 



dq k 



do-£ 0) ((^v,g<jv)ZW_i(q<Ar| • • • |<?o) 



M-l 



E *V 



k=0 



da^ iacl H^ k ,q k ) 



<lk 



dq k +i 



da 



+i 



(l.incl) , 



A k (q k+1 \q k ) A k ^(q k \ ■ ■ ■ \q Q ) 



(18) 



t>k+l, Qk+1/ 



dqk+l 



A k (q k+1 \q k ) A k _x(q k \ ■ ■ ■ \q ) 



PS 



- 1 [da^ cl \4>M,qM)A M _ 1 (q M \---\q Q ) 



which, upon parton shower action, will provide a merged sample with NLO accuracy for up to M-jet observables, LO 
accuracy for up to -/V-jet observables, including resummation at whatever accuracy is provided by the parton shower. 
Here, we have denoted the NLO corrections to inclusive cross sections as 



da( 1 ^ cl '(0„,<Z„) = daW(^,q„) + 



«» d<x(°>(^ n+1 ,g n+1 ) fl , 

dq n +i -\ 0{q n - p) 

o U<Zn+l 



(19) 



Note that the action of PS" 1 integrates to one, which clarifies once more that we are able to preserve inclusive cross 
sections. Given the extensive discussions on reclustering, dynamic scale choices and the generation of the Sudakov 
weights in the context of merging approaches so far [7H13]. we will not include a detailed algorithmic definition 
of the merging procedure here. Technical aspects will be subject to ongoing and future work concerned with the 
implementation of the procedure outlined here. 



6 Corrections below the merging scale 

Having derived the merged cross section in the presence of both tree level and one loop matrix elements, note that we 
actually have solved a NLO matching condition for each exclusive n-parton contribution above the merging scale, cf. 
eq . [T7] More precisely, 

— by preserving inclusive cross sections in the presence of only tree level matrix elements, we have matched each 
exclusive n-parton contribution to a LoopSim approximated nLO, and 

— by including exact NLO corrections, we have fulfilled this matching condition at NLO, while finally 

— by preserving inclusive cross sections in the latter case, we start to generate approximate NNLO pieces, which 
could well be the basis for NNLO matching. 

These considerations apply to contributions with emissions above the merging scale. Below the merging scale, we are 
still left with the shower approximation, and we will finally give a simple prescription of how nLO (in the case of tree 
level merging) and NLO (in the case of one loop merging) accuracy can be achieved also below the merging scale. To 
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be precise, we consider the cases of n partons above the merging scale and the first order contribution to zero or one 
emission below the merging scale. Demanding the required nLO or NLO accuracy, we obtain generalizations of the 
familiar NLO matching corrections restricted to the phase space for one emission below the merging scale, 



j ^match 



N-l 
k=0 



i (0), match 

d %\ P 



,9k)- / dq k 



da 



Ik+i - 



(0), match , 

Hp 1 



dqk 



+1 



A k (q k+1 \q k ) A k -i(q k 



+ dcr i ( |i' match ( ( / ) A'> qN)AN-i{qN I — |go) 



\lo) 



(20) 



with 



da 



(0), match / 



,q n )=daW(<p n ,q n )6(q n - 1 -p)6(p-q n )-d ( rW{ 



-1, <7ra-lJ 



p\p{ 



i-lj In) 



(21) 



Note that these corrections do not change inclusive cross sections as is the case for plain NLO matching. Within 
these contributions, the parton shower cutoff can be sent to zero provided the singly unresolved limits are reproduced 
properly; a finite shower cutoff will act similar to a phase space slicing parameter in terms of which the NLO cross 
section is reproduced. 



7 Conclusions and outlook 

We have presented an extension to multileg matrix element and parton shower merging, which preserves inclusive 
cross sections at the level of the available accuracy, particularly at tree and one-loop level. This constraint is of utmost 
importance particularly for the latter case, as NLO corrections to shower splitting kernels are so far out of reach. 
This lack of shower accuracy manifests itself in terms of potentially large logarithmic contributions which are of the 
same order of magnitude as the NLO corrections tackled in recent approaches of combining fixed order corrections 
and parton shower resummation, thus spoiling NLO accuracy for lower jet multiplicities. 

The formalism used to derive the modified algorithm is general enough to study the inclusion of even higher 
order corrections by successively replacing virtual contributions as approximated by the shower through their exact 
counterpart. The ingredients at for tree level and one loop merging, respectively, are given by fully differential LO 
calculations and inclusive NLO corrections differential in the respective Born variables. An implementation of the 
algorithm presented is subject to ongoing and future work. 



Note added 

During completion of this work a very similar approach to the problem by Lonnblad and Prestel came to the attention 
of the author, first results of which have been presented for the tree level case in [38] and for the onedoop case in [39] . 
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